ARCHIMEDEAN-TYPE FORCE IN A COSMIC DARK FLUID: 
III. BIG RIP, LITTLE RIP AND CYCLIC SOLUTIONS 
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We analyze late-time evolution of the Universe in the framework of the self-consistent model, in 
which the dark matter is influenced by the Archimedean-type force proportional to the four-gradient 
of the dark energy pressure. The dark energy is considered as a fluid with the equation of state 
of the relaxation type, which takes into account a retardation of the dark energy response to the 
Universe accelerated expansion. The dark matter is guided by the Archimedean-type force, which 
redistributes the total energy of the dark fluid between two its constituents, dark energy and dark 
matter, in the course of the Universe accelerated expansion. We focus on the constraints for the 
dark energy relaxation time parameter, for the dark energy equation of state parameter, and for 
the Archimedean-type coupling constants, which guarantee the Big Rip avoidance. In particular, 
we show that the Archimedean-type coupling protects the Universe from the Big Rip scenario with 
asymptotically infinite negative dark energy pressure, and that the Little Rip is the fate of the 
Universe with the Archimedean-type interaction inside the dark fluid. 
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I. INTRODUCTION 

The present time accelerated expansion of the 
Universe discovered in the observations of the Super- 
novae la [H-Q revived the discussion about the fate 
of our Universe. A psychologically jeopardy variant 
of the Universe future is that the permanently in- 
creasing scale factor a{t) and the Hubble function 
H{t) can reach infinite values during a finite time 
interval. The life-time of such Universe is finite, and 
the catastrophe of this type produces cosmic inertial 
(tidal) forces, which destroy the bounds in all physi- 
cal systems. Various aspects of the singular scenaria 
of the future stage of the Universe evolution were 
discussed during last two decades (see, e.g., the re- 
view Q). The specific term Big Rip (or Doomsday) 
entered the scientific lexicon after publication of the 
paper Q ; nowadays this term indicates a new trend 
in theoretical cosmology (see, e.g., It seems 

to be uncomfortable for physicist to think that the 
Big Rip is the fate of our Universe, probably, it is 
a reason that many authors consider the models, in 
which the Big Rip can be avoided. In particular, the 
Big Rip can be avoided in the F{R), /(r)-gravity 
models and their modifications, in the models with 
Chaplygin gas, in the model for the dark energy 
with various effective time-dependent equations of 
state (see, e.g., [11]- [13]). As a more optimistic vari- 
ant of the Universe behavior one could consider the 
case, when the Universe life-time is infinite, and the 
Hubble function tends asymptotically to a constant, 
H Hao- The asymptotic regime of this type ap- 
pears, in particular, in the ACDM model, which 



converts at t — > oo into the de Sitter model with 
Hoo = \J^ (A is the cosmological constant). More 
general case, when H{t) — > H^o— const, but is 
not necessary equal to \J~^, is indicated as Pseudo 

Rip in [28| . Various intermediate scenaria with the 
infinite Universe's life-time, in which the scale factor, 
the Hubble function and (probably) its time deriva- 
tive tend to infinity, belong to the class indicated by 
the term Little Rip (see, e.g., [2^-[35 for references). 



A. On the classification of the models of the 
late-time Universe behavior 

In order to classify the models of three types men- 
tioned above in more detail, we use the terminology 
(see, e.g., [lH) based on the asymptotic properties of 
the scale factor a(t), Hubble function H{t)=^, and 
its time derivative H, which are the basic quanti- 
ties in the isotropic spatially homogeneous cosmolog- 
ical models of the Friedmann-Lemaitre-Robertson- 
Walker (FLRW)-type with the metric 



dt'^-a'^(t)[{dx^f+{dx^f+{dx^)^ 



(1) 



(Here and below we use the units with c—1). In fact 
this classification is originated from the idea that 
two bonded particles are influenced by the inertial 
force proportional to the quantity 



a ■ o 



(2) 



where — g is the acceleration parameter. The inter- 
nal structure of the physical system is assumed to 
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be destroyed, when this incrtial force exceeds the 
internal (e.g., intermolecular) forces. We think that 
such force can also be indicated as a tidal force, if 
we take into account the following motives. Let us 
start with the well-known equation of the world-line 
deviation (see, e.g., [35j ) 

^ = -R'kimU'u-^n' + r , (3) 

where U^=^^ is the velocity four-vector, is 
the deviation four-vector, t is the world-line param- 
eter and A is the parameter describing the world- 
line family. The second term in the right-hand side 
of ^ describes the contribution of the force of the 
non-gravitational origin (its structure for the electro- 
magnetic interactions is described, e.g., in [s^). The 
first term in the right-hand side of ^ is interpreted 
as the tidal (curvature induced) force proportional 
to the convolution of the Riemann tensor R''^.i^ with 
the deviation four- vector. Taking into account that 
in the FLRW model the velocity four- vector is of the 
form U^=Sq, we can choose arbitrary direction, say 
Ox^ , and calculate the corresponding tidal force com- 
ponent 5R^=-i?^QjQn^ The result is ^^=^n^, thus 
the interpretation in terms of tidal force is equiva- 
lent to the one in terms of incrtial force. This force 
tends to infinity, when H or H (or both of them) 
tend to infinity. The standard classification uses the 
following terminology: 

1. Power-law inflation, if H^j, H{t ^ oo) ^ 0; 

2. ACDM model H^Hp, = where A is the 
cosmological constant; 

3. Pseudo Rip, if H{t — > oo) — > Hryo < oo; 

4. Little Rip, when H{t -> oo) — > oo, e.g., as ft-*e^'; 

5. The future singularities at finite time t^, when 
H{t ^> ts) -> 00. 

The fifth subclass can be specify, for instance, as 
follows (see [Tl|): 

5.1. Big Rip (Type I), if a ^ oo, e.g., as a — (t-\^y ' 

5.2. Sudden singularity (Type II), if the pressure 
n tends to infinity, but a{t) and energy density p 
remain finite; 

5.3. Singularity of the Type III, if 11 ^ oo and 

p — oo, but a{t) is finite; 

5.4. Singularity of the Type IV, if only iJ — >• oo. 

B. The goal of this paper 

This paper is considered to be the third (final) 
one describing the fundamentals of the model of the 
Archimedean-type interaction between the dark en- 
ergy (DE) and dark matter (DM). In [s^] we estab- 
lished the self-consistent model in which the dark 



matter is influenced by the force proportional to the 
four-gradient of the dark energy pressure. In (38| we 
classified the corresponding scenaria of the Universe 
evolution with respect to a number of the transition 
points, in which the acceleration parameter —q{t) 
takes zero values, and thus the epochs of decelerated 
expansion is replaced by the epoch of the accelerated 
expansion. In [39j we considered the application of 
the model to the problem of light propagation with 
non-minimal coupling in a cosmic dark fluid with 
an Archimedean-type interaction between the DE 
and DM, and described the so-called unlighted cos- 
mological epochs, for which the effective refraction 
index n{t), the phase and group velocities, Vph{t) 
and Vgr (t) , respectively, were imaginary functions of 
time. 

In (37l - l39| we focused on the solutions with asymp- 
totically finite Hubble function and finite accelera- 
tion parameter, and claimed that the case of infinite 
H and q is the topic of special discussion. Now we 
systematically consider the last case in terms of Big 
Rip and Little Rip solutions. Why this consider- 
ation might be interesting for readers? There are 
at least three motives to complete our investigation. 
The first motif is the following. The dark energy is 
considered in [33, HI] as a fluid with the simplest 
rheological property: the equation of state is of the 
relaxation type, i.e., it takes into account a retar- 
dation of the dark energy response, and includes an 
extra parameter, ^. The first question arises: is it 
possible to avoid the Big Rip by the appropriate 
choice of the relaxation parameter ^? The second 
motif is connected with the effective constant of the 
Archimedean- type coupling, i/* . When this coupling 
is absent, the energy density and pressure of the dark 
matter decrease in the course of the Universe expan- 
sion, so that the dark matter does not play any role 
in the asymptotic regime. However, when the pres- 
sure of the dark energy, being negative, becomes in- 
finite in the course of accelerated expansion of the 
Universe, the Archimedean-type coupling leads to 
an effective heating of the dark matter, so that the 
DE and DM contributions happen to be of the same 
order. The second question arises: for what val- 
ues of the coupling parameter v^, the Archimedean- 
type interaction is able to protect the Universe from 
the Big Rip singularity? The third motif is that 
the pressure of the dark energy, H, is the key el- 
ement in our model, which guides the behavior of 
both state functions of the dark matter, the energy- 
density E and pressure P (via the Archimedean-type 
force), the energy-density p of the dark energy (via 
the equation of state), thus predetermining the be- 
havior of the Hubble function H and its derivative 
H (via the Einstein equations). This means that 
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in the model under consideration one can classify 
the asymptotic types of the Universe behavior us- 
ing only one element, just the function n(t). This 
function, n(t), satisfies the key equation, which is 
the nonlinear differential equation of the second or- 
der and includes six guiding parameters: ^, a and po 
coming from the three-parameter equation of state 
for the dark energy, the effective coupling constant 
of the Archimedean-type interaction z/* , the starting 
energy, i?*, and temperature, T*, of the dark matter. 
Thus, (keeping in mind the initial data n(<o) and 
n(to)) we have a possibility to classify all the types 
of asymptotic behavior of the model using eight pa- 
rameters. The question arises: in what domains of 
the 8-dimensional effective space of the guiding pa- 
rameters the Big Rip is reliably avoided? 

In order to answer these questions we organized 
this paper as follows. In Section II we remind 
briefly the master equations for the DE and DM 
evolution in the framework of the Archimedean-type 
model (Subsections IIA-IIC). In Subsection IID we 
describe our classification scheme, which is based 
on eight parameters and includes three bifurcation 
points: i^*=0 {v^, is the Archimedean-type coupling 
constant); ^=0 (^ is the DE relaxation parameter); 
0"=— 1, the phantom-crossing point. In Section III 
we focus on the analysis of the models indicated as 
models with asymptotic DE domination, for which 
the contribution of the DM to the Universe expan- 
sion is asymptotically negligible. In Subsection IIIA 
we describe the Big Rip, Pseudo Rip, Cyclic solu- 
tions with v^:=Q and ^=0 in terms of two guiding 
parameters ct, poi and the initial value of the DE 
pressure. In Subsection IIIB we focus on the mod- 
els with i/*=0 and ^ ^ 0, thus analyzing the role 
of the DE relaxation parameter in various scenaria 
of the Big Rip and Little Rip formation (in terms 
of three guiding parameters ^, a, po-, and two ini- 
tial data for the DE). In Section IV we consider 
the models with 7^ and ^=0 (in terms of three 
guiding parameters , a, po and the initial value of 
the DE pressure), and show that the Archimedean- 
type coupling of the DM with the barotropic DE 
avoids the Big Rip regimes in the Universe expan- 
sion. In Section V we discuss the general case i^* 7^ 
and ^ 7^ qualitatively and numerically, and show 
that the Archimedean-type coupling of the DM with 
the DE of the rheological type converts the Big Rip 
regimes into the Little Rip, Pseudo Rip and Cyclic 
regimes of the Universe expansion. 



II. MASTER EQUATIONS 

Let us remind briefly the master equations of the 
model of the Archimedean-type coupling between 
the dark energy and dark matter. These equations 



are derived in [37| and now we use them to analyze 
the late-time behavior of the Universe in the frame- 
work of this model. 



A. Equations for gravity field 

The FLRW-type isotropic spatially homogeneous 
cosmological model with two interacting con- 
stituents: the dark energy and dark matter are 
known to be described by two nontrivial Einstein 
equations 



ii = -AttGIp + E + II + P] 



SttG 



(P + E)- 



(4) 



(5) 



Here the dot denotes the derivative with respect to 
time, the functions p{t) and U{t) are the energy- 
density and pressure of the dark energy, respectively, 
E{t) and P{t) describe the corresponding state func- 
tions of the dark matter. The cosmological con- 
stant A is considered to be incorporated into the 
DE energy density p and DE pressure 11. As usual, 
H{t) = ^ is the Hubble function. The conservation 
law as the compatibility equation for the set (IH)-® 
has the form 



P + E + 'SH{p + E + II + P) = 0. 



B. Dark matter description 



(6) 



The energy density and pressure of the one- 
component dark matter can be effectively presented, 
respectively, by the integrals 

p roo 

E{x)^^ q^dq^l+q^F{x) e^^'V^^, (7) 



P{x) 



^E^ r F{x)qHq 



3x3 ^1+q^Fix) 



(8) 



The auxiliary function F{x) of the dimensionless 
variable x = -^7^ 

a(to) 



F{x) = — exp {2^., [H(1)-H(.T)]} , (9) 
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and convenient parameters 



(10) 



contain the effective Archimedean-type couphng 
constant j/*, an effective number density iV*, the 
mass m* and temperature T* of the leading sort of 
the dark matter particles. The term 

/•OO 

Ks{\^,) = / (i^icoshsz • exp [— A* coshz] , (11) 
Jo 

is the modified Bessel function, fc(B-) is the Bohz- 



mann constant. 



C. Dark energy dynamics 



In these terms the balance equation ^ and the con- 
stitutive equation (jl2l) give the key equation for the 
pressure of the dark energy n(x) 

^x^U"{x)+xU'{x) (4^+(7) +3{l+a)U+3po^J{x) . 

(16) 

The prime denotes the derivative with respect to x. 
The source J^{x) is defined as follows: 



^iWFix) 



The quantity J^{x) vanishes, when the Archimedean- 
type coupling parameter vanishes, i.e., z/* = 0. We 
deal with the differential equation of the second- 
order linear in the derivatives but nonlinear in the 
unknown function Il{x). 



To describe the dark energy fluid we use the linear 
three-parameter equation of state 



p{t) = po + all 



Hit) 



n. 



(12) 



When a=0 and ^^=0, Eq. introduces the model 
in which the dark energy relates to the cosmological 
constant A. When po=0 and ^=0, Eq. (fT2|) gives the 
well-known linear relation 11=11} p with w = The 
retardation of response in the dark energy evolution 
is taken into account by inserting the term contain- 
ing the first derivative of the pressure 11. An equiv- 
alent scheme is widely used in the extended ther- 
modynamics and rheology (see, e.g., [405), in which 
the extended constitutive equation for the thermo- 
dynamically coupled variables X and Y has the form 



tX + X = w Y . 



(13) 



Here t is a relaxation time, a new coupling param- 
eter of the model. In the cosmological context r is 
generally the function of time, T{t). We assume that 
^(^)~CT^(t)' ^^^^ relaxation time can be mea- 
sured in natural cosmological scale [4l| . Our ansatz 
here is that the dimensionless parameter ^ is con- 
stant. 

When the quantities p{t) and Il{t) depend on time 
through the scale factor a{t) only, i.e., p=p{a{t)), 
n = n(a(i)), the so-called x- representation is con- 
venient, which is based on the following relations 



d 
dt 



t-tn = 



xH{x) 



dx ' 



i(to) 



dx 
xH{x) 



(14) 



(15) 



D. Classification scheme based on the analysis 
of the DE pressure 11 

The function 11 plays the key role in the analysis 
of the whole model behavior. Indeed, when n(x) is 
found from the key equation (|16p . we can reconstruct 
the DE energy-density p using the formula 



p[x) = po + crn(a;) + S^xU' [x) 



(18) 



then calculate the DM state functions E{x) and 
P{x) using (O and dS]), then find H{x) from the 
Einstein equation ([5]), and finally, reconstruct the 
scale factor as the function of time using (fT5|) . 

We are interested to discuss the late-time period 
of evolution, i.e., the period, when t ^ too- We con- 
sider both cases: too=ts (future finite time singu- 
larity) and too=oo. Clearly, there are three possible 
types of the asymptotic behavior of the dark fluid 
composed of the dark energy and dark matter. 
1. The first type refers to the asymptote n(i — > 
^oo) — ^ Iiao=const. This case includes the submodel 
with noo=0). Since the Archimedean- type force is 
proportional to the four-gradient of the dark energy 
pressure 11, then asymptotically the dark matter de- 
couples from the dark energy, its energy-density E 
decreases as and we obtain the case indicated as 
an asymptotic DE domination. 

2. The second type of behavior refers to the case 
n(t too) ^ +00. Now the DM becomes frozen, 
since the function F(x) ^ tends to zero exponen- 
tially, the DM decouples from the DE, and again we 
deal with the model of the DE domination. 

3. The third and the most interesting case, refers 
to the asymptotic behavior H{t — > ioo) ^ — oo, for 
which the function F{x) tends to infinity, the DM 
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becomes effectively ultrarelativistic and tlius plays 
an active role in the energy redistribution process. 

The model as a whole contains four effective guid- 
ing parameters a and po describing the DE equa- 
tion of state, and j^*, the effective Archimedean- 
type coupling constant) and four initial parameters 
(n(to) and n(io)j initial data for the DE pressure, 
and i?*, r* initial data for the DM). The classifica- 
tion scheme includes three bifurcation points: the 
first, i/*=0 indicates whether the Archimedean-type 
coupling is switched on or not; the second, ^=0 in- 
dicates whether the DE possesses the simplest rheo- 
logical property or not, the third, (t=— 1, relates to 
the phantom-crossing point. 

Remark on the Big Rip symptom: 

when the integral in the right-hand side of (|15p con- 
verges on the upper limit, the infinite value of the 
scale factor a{t) can be reached at finite time t(s)j 
thus we deal with the Big Rip. 



III. BIG RIP SOLUTIONS IN THE MODELS 
WITH ASYMPTOTIC DARK ENERGY 
DOMINATION 

A. Big Rip solutions at ^=0 and !/*=0 

This first model relates to the case, when the dark 
energy does not possess rheological properties, and 
there is no Archimedean-type coupling between DE 
and DM. The contribution of the DM is asymptot- 
ically vanishing, and in this sense we deal with the 
case of the DE domination. The main results for 
this model are well-known and our goal is to recover 
them as the limiting case of our model in the corre- 
sponding terminology. 



1. Power-law expansion with po=0 

Let us recover, first, the well-known results pre- 
senting the Big Rip solutions, when ^=0, po=0 and 
^^=0. The key equation takes the form 



and the solution is 



n(x) = n(i) X 



(19) 



(20) 



It is convenient to introduce the new dimensionless 
parameter 



'3(l+a) ' 



(21) 



since it plays an important role below. For the DE 
energy-density p(x) we immediately obtain 



so that 



p{x) = crn(l) x^' 



H^{x) = -87rGcrn(l) 



(22) 



(23) 



Clearly, the real solution for the Hubble function ex- 
ists, when the product crn(l) is positive. Using the 
remark on the Big Rip symptom one can conclude 
that the integral in (jl5p converges, when a is pos- 
itive, i.e. — 1 < (T < 0. This simplest submodel 
allows us to verify this fact directly by computing 
the scale factor analytically. Indeed, according to 
(flSl) the scale factor has the form 



a{to) 



l-:^V67rGan(l)(t-io) 
da 



(24) 



When — 1 < cr < or in other words, when w=— < 
— 1, one obtains that the parameter a is positive. 
This means that the scale factor tends to infinity at 
some finite time value i=t*, given by 



3a 



^67rGcrn(r 



> to. 



(25) 



When cr < 0, the solution is real, when the initial 
value of the DE pressure is negative 11(1) < 0. Thus 
we deal with the well-known Big Rip solution, which 
is characterized by the following asymptotes a.t t ^ 
t^,: n — — oo, p — )■ -l-oo, a — )■ oo, H ^ oo and 
ij — )■ oo. Let us add that when the parameter a is 
negative and belongs to the interval — 1 < cr < 0, 
the parameter a is positive (see ([IT])). There are 
two special cases. When cr — > — 1, the parameter 
a becomes infinite, a — > +oo; when cr — >■ 0, the 
parameter a tends to zero also. Finally, when a is 
negative, i.e., a > or cr < —1, the solution ([M)) 
describes power-law expansion of the Universe. The 
Little Rip solutions do not appear in this submodel. 



2. Solutions with po ^ 

When po ^ 0, the solution of the key equation 

CTxn'(x)+3(l+o-)n+3po=0 (26) 
has the form 



n{x)-- 



1+0- " 



n(i)- 



Po 

1+a 



(27) 
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The square of the Hubble function reads now 



H\x)-- 



'3(l+a) 



where 



A = (J 



(l+a)n(l) 



(28) 



(29) 



The signs of the parameters A, pa, a and a provide 
important details of the classification. In order to 
classify the models we, first, distinguish three (prin- 
cipal) cases and then show that all other cases can 
be reduced to these three ones. 

(1) The model with A>0,-l<a<0, po>0- 
Let us consider, first, the case when — 1 < <t < and 
A > 0. The corresponding initial value 11(1) should 
be negative, and the inequality 11(1) < —-^^ has to 
be satisfied. The solution for the scale factor has the 
form 



a{to) 
where 



cosh fi{t— to)— A sinh \p\{t—to 



p : 



■^6TTGpoil+a) 



(30) 



(31) 



There exists a moment t—t^, for which a(t*)=oo. At 
t t^ the plot oia{t) has a vertical asymptote. This 
critical moment of time, t^, can be found from the 
equation 



tanh \p\{t^,—tQ) 



1 



VTTa 



and has the form 



1 , fvnTA+i\ ^ 



(32) 



(33) 



Again we deal with the Big Rip solution, character- 
ized by the asymptotes IT — > — oo, p — > -t-oo, a — > oo, 
— oo, ij — ^ oo, but now, in the case po 0, 
we can indicate the Big Rip singularity as that of 
hyperbolic type in contrast to the one of the power- 
law type, presented in the previous subsection (with 
Po=0). 

Remark 1: 

Two submodels with negative a, namely, A> cr > 
0, Po > and A> a < —1, po < 0, give formally 
the same scale factor However, the essential 

difference is that at a < the Hubble functio n H 
(1^ tends to the constant H -^00 = ^^^^, and 
we deal with the Pseudo Rip instead of the Big Rip. 



(2) The model with A < 0, -I < a < Q, po > Q. 
When — 1 < (T < and ^ < 0, one obtains from 
that there is a critical value x—x^ 



(34) 



for which the function H^{x) 
iirGpa 



H\xy 



3(1+0-) 



1- — 



(35) 



takes zero value. Let us assume that x^ > 1, i.e., 
according to dMl) < 1 and -I < A < Q. It 
is possible when --pq^ < n(l) < Since the 
model solutions for H{x) obtained from ([35]) can not 
be prolonged for x > x^,, x starts to decrease after 
the moment i=t(max)- The corresponding solution 
for the scale factor has now the form 



a{to) 



cosh /i(^o-^ (max)) 
C0Sh//(i-t(max)) 



(36) 



where 



W)-^o + ^log(^i±^3j>t„. (37) 

Clearly, the plot of this function is symmetric in ref- 
erence to the moment i(niax) ] the formula (j36p gives 
the maximal value (p4)) a{tt)-=a{to)\A\~°' at this mo- 
ment. At i ^ 00 the scale factor a{t) tends to zero 
as 



a{t) cx expj-iJooi} , H^. 



3(1+ct) 



(38) 



The Hubble function is described by the finite func- 
tion 



H{t) = -iJootanh [|Ml(i-t(max))] 



(39) 



which changes the sign at t=t 



(max) ■ 



This means that 



the Universe expansion turns into a collapse at this 
moment. The acceleration parameter 



'9(i) = ^^ = l-^sinh V(i-^(max)) (40) 



takes infinite value at t=i(niax) i since the Hubble 
function in the denominator vanishes at this mo- 
ment. Clearly, the Universe passes two eras with 
—q>0 (accelerated expansion and accelerated col- 
lapse), and two eras with —q<0 (decelerated ex- 
pansion and decelerated collapse). We deal with a 
closed Universe which does not admit the Big Rip 
behavior, despite — 1 < cr < 0. 
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Remark 2: 

Two submodels with negative a, namely, ^ < 0, cr > 
0, po < and A < 0, a < —1, po > 0, give formally 
the same expression for the scale factor Since 
now a is negative, the Hubble function H (1^51) tends 

again to the constant H — !• -ff oo — \J'' t{i+'a) ' 
require < 1 to avoid the internal singular point. 
We deal now with the Pseudo Rip instead of Big 
Rip. 

(3) The model with A < 0, -I < a < 0, pq < 0. 
Now the scale factor can be expressed in terms of 
trigonometric functions 

Q^U) r 1 -2a 

cos\p\{t-ta)-^\A\-lsm\n\{t~to) 
a[to) L J 

(41) 

the parameter A has to satisfy the inequality A < 
— 1, thus, the initial value 11(1) is restricted by 
n(l) < —■'i^- Again, there exists a moment 



1 

to + f-rarctg 



1 



(42) 



for which the plot of the scale factor has the vertical 
asymptote typical for the Big Rip with 11 — — c», 
p — !■ +CX), a — > oo, H — !■ oo, H — > oo. It is con- 
venient to indicate such singularity as the Big Rip 
singularity of the trigonometric type. 

Remark 3: 

Two submodels with negative a, namely, ^ < 0, 
cr > 0, Po < and ^ < 0, cr < —1, po > 0, give 
formally the same scale factor (|4T|) . Since now the 
parameter p becomes imaginary (see (PT|) ). at a < 
one should replace the trigonometric functions by 
the hyperbolic one. Again, asymptotically we deal 
with the Pseudo Rip solution instead of Big Rip. 



B. Big Rip and Little Rip solutions at ^ 7^ 

and !/*=0 



When we consider this submodel we assume that 
the Archimedean-type coupling is switched of, but 
the dark energy possesses the simplest rheological 
properties. Since j^*=0, the asymptotic regime again 
will be characterized by the dark energy domination. 
In order to analyze this model let us introduce the 
new unknown function Z{x) as follows: 



n(x) 



-Z{x) , 



and define the initial value Z(\) as 

Po 



z(i) = n(i) 



(43) 



(44) 



The function Z{x) satisfies the Euler equation 

ix^Z"{x) + xZ'{x) (4C + 0-) + 3(1 + a)Z = 0. (45) 

The characteristic polynomial of this Euler equation 
has two roots 

-1.2 = ^ 



-{a + 2,0±V¥^W^M , (46) 



which can be real or complex depending on the val- 
ues of the parameters ^ and a. One can distinguish 
three subcases. 



1. Two different real roots [{a - 3^)^ > 12^] 

When the discriminant in (j46p is positive, one ob- 
tains 

Z{x)^x-^ [^{x^+x-^) + 



[Z'{l)+jZ{l)] ( r-T, 



(47) 



where 



7^^^, r^lv/('^-3e)2-12e. (48) 

According to the remark about the Big Rip symptom 
we can conclude that the integral in (fTSl) converges 
now, when F > 7, thus this inequality guarantees the 
Big Rip existence. In order to illustrate the existence 
of the Big Rip by exact solutions, let us assume that 
initial data satisfy the equality 



z'(i) = z(i)(r-7) = (r-7) 



n(i) 



Po 



l+a 



(49) 



Then the solutions can be simplified as follows 

Z{x)^Z{l) x^-'' , (50) 

and Z{x) becomes proportional to the required lead- 
ing order term exactly. Now the DE state functions 
are of the form 



Il{x)-- 



Po 

l+a 



-Z(l) x^-T. 



p(x) 



Po 



n{^)z{i) x 



r-7 



1+cr 

where the new parameter 

n{0 = a+ar-i) 



(51) 



(52) 



(53) 
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is introduced. In the case of the DE domination the 
square of the Hubble function 



'3(l+cr) 



i^mii^ z(i) 

Po 



(54) 

formaUy coincides with (|28l) , if we use the substitu- 
tions 

^^m^zii), i^r-,. (55) 

Pa a 

Based on the resuhs of the previous section one can 
conclude that the Big Rip takes place, if F > 7. 
There are two domains on the plane of the parame- 
ters (tO^, where this inequality holds at ^ > 0: first, 
-I < a < C < 5, second, <t < -1. The 

first domain is a natural extension of the interval 
— 1 < a < 0, which is obtained if we put ^=0 into 
the right-hand side of the inequality; let us remind 
that just for the case — 1 < u < the Big Rip ap- 
pears, when ^=0. The second interval for the Big 
Rip existence, cr<— lat^>0, has no analog at 

The parameter po distinguishes three situations 
admitting the Big Rip at ^ > 0: 

(1) Po=0. 

In both domains the singularity of a{t) is of the 
power-law type. 

(2) Po > 0. 

In the first domain, where — 1 < tr < 3^—2^/3^, 
^ < i, the singularity is of the hyperbolic type. In 
the second domain, where cr < — 1 , the singularity is 
of the trigonometric type. 

(3) Po < 0. 

In the first domain, where —1<(j< 3£,—2^/3^, 
^ < i, the singularity is of the trigonometric type. 
In the second domain, where ct < — 1, the singularity 
is of the hyperbohc type. 

More detailed constraints for the guiding param- 
eters follow from the inequalities rt{^)Z{l) > or 
ri(^)Z(l) < 0, respectively. 

(it) C < 0. 

When the parameter ^ is negative, we use the fol- 
lowing scheme of analysis. We put 



z'(i) = -z(i)(r + 7) 

and obtain that 

Z{x)=Z{l) x-'^^+^K 
The DE state functions take now the form 



n(x) 



Po 
l+a 



+Z{1) X 



-(r+7) 



(56) 



(57) 



(58) 



Po 

l+a 



(59) 



and the square of the Hubble function reads 



SirGpo 
3{l+(j) 



»(0(i+'t) 

Po 



Z{1) X 



-(r+7) 



(60) 

Again, these functions can be obtained asymptoti- 
cally at X 00, since now F < and in this limit 
one obtains that — !• 0. The sum r-l-7 should be 
negative, and the product r2(^)Z(l) should be pos- 
itive, if we try to find the analogs of the Big Rip. 
This is possible, when a > —1. Clearly, we obtain 
the power-law singularity, when po=Oj the hyper- 
bolic singularity, when po > 0, and the trigonomet- 
ric singularity, when po > 0- 



2. Double real roots [{(t-3^)^ = 12^] 

The discriminant can take zero value, only when 
f > 0. For the positive ^ there are three interesting 
cases: ^ < ^, ^ > ^ and C=|- 

The corresponding solution for Z(x) is 

Z{x) = x-^ {Z{1)+ log X [Z'{l)+jZ{l)]} , (61) 

where 7=3+^1, if a^3£_+2y/3^, and 7=3-^1, if 



a=3^—2^/3^. According to the mentioned symptom, 
the Big Rip is possible, when 7 is negative; clearly, it 
is possible only if we deal with the second solution, 
i.e., when^ < i and cr==3^— 2^/3^. These constraints 
require that — 1 < u < 0. In order to illustrate the 
behavior of the scale factor, let us assume that, first, 
po=0, second, Z'(l)=-Z(l)f and, third, Z{1) < 0. 
Then we obtain 



p{x) = 3|Z(1)| xl'^lloga;. 



l7h^^(V3e-l). 

The square of the Hubble function reads 
H^{x) =87rG|Z(l)|a;l'^llogx, 
and the integral in ([15]) 



^|7| log a; 



V4G|7Z(l)|(<-to) = erf 

happens to be reduced to the error-function 

erf [z] = —= / du e"" . 
Vi" Jo 



(62) 



(63) 



(64) 



(65) 



(66) 
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The error-function erf [z] takes finite value erf [oo] = l, 
when the upper limit in the integral goes to the in- 
finity, thus, the value x — oo can be reached at the 



moment ts 
form 



--to+ 



The scale factor has the 



a{t)=a{to) exp <j — 



d-'WAG\jZmt-to)] 



(67) 

where the symbol erf ^ stands for the function in- 
verse to the error-function. Thus, when ^ < ^ the 
Big Rip is possible, if a=3^—2^/3^. 

This case corresponds to the positive parameter 7. 
This means that Z{x) (see ((6T|) ) asymptotically van- 
ishes, the Hubble function tends to constant, and we 
deal with the Pseudo Rip. 

(Hi) C = 3- 

Double real roots appear now, when (ct— 1)^=4, i.e., 
when (7=3 or (t=— 1. The first root, a— 3, corre- 
sponds to the solution (j6ip with 7= 6, and we again 
deal with the Pseudo Rip. The second root, a=—l, 
refers to the special case with diverging quantity 
, and we have to return to the key equation . 
As we emphasized in (37j . in this special model the 
DE pressure and the DE energy density contain the 
logarithmic terms squared. In particular, for the ini- 
tial data, which satisfy the equality p(l)-|-n(l)=0, 
these DE state functions are 



n(a:)=n(l)-3po log x-^Po log^ X , 



p{x) = p{l) + -polog X. 



(68) 



(69) 



The scale factor has the form 
ait) ' ' 



a(io) 



: exp • 



'Milsinh 
9/30 



(70) 

it describes the super-exponential expansion. The 
Hubble function is monotonic 



Hit). 



87rGp(l) 



cosh 



it-to) (71) 



and increases infinitely. The function H 



iJ(0=47rGV2/9o/5(l)sinh y/UnGpo it-to) (72) 

00. The acceleration 



also tends to infinity at t 
parameter 



-qit) = 1 + . 



9po sinh [V127rG>o (i-^o)] 



2p(l) cosh^ [V127rG>o (i-io)] 



(73) 



starts with — g(fo)=l, reaches the maximum 

=1 



9po 



^-TTT at t=to+ (^+^) and tends 

y/l2nGpo 



asymptotically to — g(cx))=l. This is the explicit ex- 
ample of the Little Rip. 



3. Complex roots /(cr - 3^)^ < 12^7 

Complex roots appear at positive ^. The solution 
is quasiperiodic 

Zix) = X'"' {Zil)cosi(3\ogx) + 
[Z'il)+jZil)] 



where 



7 ^ 



3^ 



■sin(/31oga:) L (74) 



1 



-ia- 3^)2 . (75) 



2C ' - 2^ 

In this section it is convenient to use the variable 
T= log X and to present the square of the Hubble 
function as follows 

BttG/Oo 



H\r)= 



3(1+0-) 



{ 1+e"'^^ [A cos pT+B sin ^t] } 



(76) 
(77) 



where the parameters 

A = (TZil)+iZ'il). 

B=-^i<T- [^'(1)+7^(1)] - iPZil) (78) 

depend on the initial data ^(1) and Z'il)\ for the 
illustration one can choose them so that S=0. Again 
there are three interesting cases. 

(i) H'^ir) > for arbitrary t. 

For instance, when 7 > 0, B—0 and \A\ < 1 we 
obtain that H^ir) > 0. Then the Hubble function 

Hir) is real and tends to Hco 

we deal with the Pseudo Rip. 



3(1+<t) '^'^ ^ 



00; 



(a) H^ir) > 0. 

To illustrate the situation, for which is nonneg- 
ative and has a number of zeros, let us choose, first, 
that 7=0, second, that A=l, B=0, providing 



n , , 167rG/9n o 1 



(79) 



The Hubble function is periodic with maximal value 



' (max) 



■ The corresponding scale factor 



a(t)=a(to) exp < — arcsin 



(80) 
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takes finite values satisfying the inequalities 

ait) 



e f> < 



a{to) 



< ef> . 



(81) 



Thus, there are neither Big Rip, nor Little Rip, de- 
spite -1 < cr < < 5). 

(in) changes the sign. 

When 7 < 0, the quasiperiodic function H'^{t) given 
by ()76p inevitably reaches zero value at some mo- 
ment T=T^,. Near this point the function H'^{t) can 
be presented as 

~ /i(r,-T) + icj(r,-T)2 + ... (82) 

Depending on the guiding parameters a, ^, poj a-nd 
on the initial data Z{1) and Z'{\), the quantity h 
can vanish or be positive. 
{l)h^Q.^ . 

In the vicinity of the moment one obtains that 
r»-r=|(t-r)2, and thus a{t) oc exp [-^i{t-t*f] . 
(2) /i=0. 

The decomposition (|82|) starts with (r*— t)-^, we ob- 

ri_e-y¥(*-*-)" . 

In both cases the Universe has finite size and its 
life-time is finite. 



tain that a(t) oc exp 



4- Resume 

To summarize the results of the analysis let us in- 
dicate the domains on the plane of the parameters 
^ and a, for which the Big Rip scenaria can be re- 
alized: these domains are displayed on Fig.l. First, 
let us focus on the case ^=0. According to the re- 
sults of Subsection IIA, the Big Rip can be realized 
on the interval — 1 < cr < of the vertical line 
It is the "classical" Big Rip domain, which is char- 
acterized by the condition p-\-\i=^^p < 0. When 

< ^ < i this "classical" Big Rip zone contracts 
along the line a to the interval — 1 < cr < 3^— 2-\/3f, 
and disappears at all when ^ > ;5 (see the domain 

1 on Fig.l). In this sense, the simplest rheological 
property of the DE, namely, the retardation of the 
DE response to the Universe expansion, can avoid 
the Big Rip, which was the fate of the Universe at 
— 1 < CT < 0, if the relaxation parameter ^ exceeds 
the critical value £,=\- In other words, when the DE 
relaxation time T{t) = is bigger than 377=^, 

where Q = "^kU^ is the Universe expansion scalar, 
the regime of the Big Rip can not be supported by 
the DE with such relaxation time. If ^ is negative, 
the interval — 1 < cr < again relates to the Big Rip 



scenario of the Universe expansion (see the domain 

II on Fig.l). The choice of the parameter po can not 
avoid the Big Rip scenario, but it predetermined the 
type of the Big Rip: the power-law type, hyperbolic 
or trigonometric ones. 

When cr > 0, there was no classical Big Rip at 
^=0. The same fact can be indicated at ^ > 0. 
However, if the parameter ^ is negative, the Big Rip 
scenario is possible for cr > 0, since instead of re- 
tardation the DE displays the acceleration of the 
response to the Universe expansion (see the domain 

III on Fig.l). 

When cr < — 1, there was no "classical" Big Rip 
at ^=0. However, the Big Rip becomes possible at 
cr < — 1 for positive relaxation parameter (see the 
domain IV on Fig.l). Mathematically, this new do- 
main for the Big Rip scenaria appears, since the 
characteristic equation is quadratic, and the pres- 
ence of the second root in (^Sl) extends the possibil- 
ities of the Big Rip scenaria. Indeed, when ^ 7^ 0, 
but ^ — 0, the roots of the mentioned characteristic 
equation can be estimated as 



Sl 



,(l + a) _ 1 



S2 



cr 
1 



(83) 



The first root does not depend on and gives us 
the "classical" Big Rip, when -1 < cr < 0. The 
second root is positive and thus describes the "non- 
classical" Big Rip, when either cr < and ^ > 0, or 
cr > and f < 0. This explains the appearance of 
two new domains IV and V on Fig.l. For small val- 
ues of ^ the second characteristic number S2 — ^ — ^ 
is respectively big; from the physical point of view, 
we can speak about instability of the DE response 
to the Universe expansion near the bifurcation line 
?=0. 



IV. BIG RIP AVOIDANCE AT v, / 0, C=0 

In this model we suggest that the dark energy is 
described by the equation of state p — po+cH, but 
now the Archimedean-type coupling is switched on. 
In other words, there are two coupled energy reser- 
voirs in the Universe: the dark energy and dark mat- 
ter, and an effective energy redistribution between 
them is possible in the course of the Universe accel- 
erated expansion. 



A. Asymptotic behavior at 11 — > —00 

Let the moment t^ (or equivalently, x—1) be cho- 
sen so that the value H(l) is large and negative. The 
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FIG. 1. The domains on the plane of the parameters 
^ and (T, for which the Universe evolution follows the 
Big Rip scenario in the model of the asymptotic dark 
energy domination. The horizontal stripe — 1 < cr < 
is divided into three domains. There are two domains of 
the Big Rip in this stripe: the first domain is bounded 
by the part of parabola cr^S^— 2-/3^ (0 < C < |)i by the 
vertical line ^=0 and horizontal line a—~l; the second 
domain is situated at ^ < 0. Two rectangular sectors 
5 > 0, (7 < — 1 and ^ < 0, (T > also refer to the case of 
Big Rip. 



1. New dynamic variables 

Let US introduce new variables X, Y and param- 
eter a*. For the case a > —1 the transformations 
have the form: 



n 



Po 



l+a 



(87) 



3(l+a) 



When (T < —1, we replace formally l+a by ll+crj 
and X by ~X. For the case a——l the mentioned 
transformations are not appropriate and we consider 
this submodel in the special section. In terms of new 
variables the key equation can be transformed into 
the linear equation for the function X(Y): 



Y^=AaX + e' 
dV 



(89) 



The energy density of the DE and DM in these new 
terms can be written as 



1 + cr 



cr 

— Y. 



E 



Aav^X 



(90) 



thus, the square of the Hubble function takes the 
form 



8nG 



Po 
l + a 



Y 



1 



AaX 



(91) 



term J (see (fT7|) ) takes the form 



e-'-nd) n'(x) e-^'-nix) ^ 



(84) 



where we use the new parameter 



q dqe 



A. 







A? 



(85) 

Here and below the symbols with asterisk: /*, i?*, 
^/* , relate to the parameters of the leading DM com- 
ponent. The key equation (|16l) transforms into non- 
linear differential equation of the first order 



xH' 



E^v*!* [n(i)-n]_ 



-3 [po+{l+a)Il] 



(86) 

and its structure prompts us the following transfor- 
mations. 



and the rate of the Hubble function evolution H is 
described by the formula 



H = 



AirGil+a) 



X 



(92) 



The function a{t) can be extracted from the equation 



1287rG 



Y{t) 



dY 



4a- 



X(Y) 



(93) 

Here we have to put X{Y) as the solution of 
and to use the inverse function 



Y{t) = Y[X{t)] = Y 



ait) 



^a(to)a* 

on the upper limit of the integral. 



(94) 
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2. Qualitative analysis: divergence of the basic integral 



3. Search for the scale factor 



The solution to the equation (|89p is 

4a 



X{Y)^X{1) 



Y 

F(T) 



Y 



4a 



Y(l) 



du 

,1+4q' 



(95) 

Here the constants X{1) and 1^(1) are defined using 
the initial data as follows: 



X{1) 



1 



Y{1) = -V, 



n(i) 



l + <7 



(96) 



Whatever the parameter a is, when 11 — > — oo (i.e., 
Y — >■ +(X)) the leading order terms in the asymptotic 
decomposition of the function X{Y) 



X{Y) 



1 



l+4a (l+4a)(2+4a) 

1 + ^^ + ^ TTT. - + ■ 



Y 



come from the generalized integral exponential 



(97) 



Y(X) 



du 

,1+4q' 



(98) 




which appears in the right-hand side of (|95l). Thus, 
when X — > yc^ , the integrand in behaves at 
y — )■ oo as 



(99) 



This asymptotic estimation is valid when cr < 3, i.e., 
in both interesting cases: — 1 < ti < and a < —1. 
Big Rip can exist, when the integral in the right- 
hand side of (|93|) converges at F ^ oo. In our case 
the integral diverges at F — >■ oo as VY, and Big Rip 
can not be realized in the scenario with ^=0 and 
i^* 7^ 0. In the asymptotic limit Y oo the formula 
2]) for the yields 



BttG 



1 + C7 



(3-<t) 



Y 



(100) 



In other words, when (t < 3 the Hubble function 
infinitely grows as -> \[^^^^ ■ As for H, ^ 
yields 



H 



'4 



i.e., the rate of the Hubble function growth tends to 
the positive constant Hoo—^^^^^~—^- 



Keeping in mind (l93t and (1991) we can reconstruct 
now the function a(t) using the function 



Yit) 



?7rG(3 - 



3i^* 



(102) 



for the case < a < oo. When F — oo, the behav- 
ior of the function X{Y) is predetermined by the 
second term in (j95p . since the generalized integral 
exponent gives the exponential leading order 
term X (x ye^ (see (|97l) ). The inverse function 
Y{X) is of the logarithmic type, which can be pre- 
sented by the following iteration procedure 



Y log XY log [X \og[X log[X...]]] , (103) 
zero-order estimations yielding 

Y ^\ogXY{l). (104) 

The scale factor 



a(to)a* f27rG(3-cr) 
y4(i) [ 3i^* 



(i-io)' 



(105) 



is described by the anti-Gaussian function, which 
was obtained in 37| as the exact solution of master 
equations of the Archimedean-type model. Since the 
corresponding Hubble function 



H{t) 



47rG(3-cr) 

3^A^ 



(t-io) 



(106) 



grows linearly with time, the obtained solution can 
be classified as the Little Rip with H ^ oo and 
H — > const. 

To complete the analysis let us consider two spe- 
cial cases. 

(i) Special case a=—l 

In the special case: a=oo, or equivalently, a=—l, we 
should consider the solutions in more details. The 
key equation 



,z..[n(i)-n]^-^ 



=3po 



can be now transformed into 
dX 4 



X 



using the modified replacements 



X EE 



(107) 



(108) 



(109) 
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f = -z/, [n - n(i)] , f(i) = o. 

The solution to (ITU5| 



(110) 



X{Y) 



1(1)- 



1+ 



1 + 



3po i/. 



ill!) 

shows that X{Y) takes zero value at some Y=Y^, 
and then changes the sign. Since X is a positively 
defined function, we can conclude that the regime 
a — >■ oo with n — > —oo can not be realized, when 
(T=— 1. In other words, the phantom-crossing value 
(T=— 1, which corresponds to a^oo is the critical 
value: when — 1 < cr < 0, the Little Rip regime 
with a{t) — > oo and Il{t) — oo is possible; when 
(T=— 1, the scale factor can not reach infinite value. 

(ii) Special case (7=3 

The special case: a——j, or equivalently, a=3, re- 
lates to the ultrarelativistic DE, when pq=0, since 
p=3n. The functions X, and H take now the 
form 



2ttG 



12KY 



Po 



H = 



iy,iK + eY) 
IQttGK 



[K + e^) 

where the new constant K is introduced 
K = X{\)Y{\)-e^(^K 

When y ^ oo, H - 



(112) 
(113) 
(114) 

(115) 



^^,ij^O, thus, we deal 
with the de Sitter asymptote. 



BIG RIP AVOIDANCE IN GENERAL 
CASE: C / AND ^ Q 



In terms of the dynamic variables ([87]) the key 
equation for the DE pressure takes the form 



16^ 

3(H-C7) 



X'^Y"+Y' 



4(.+7e) 



-y = 0. 



3(l+a) 

(116) 

At ^=0 it coincides with as it should be. The 
function H and the square of the Hubble function 
can be presented as follows 



H 



4nG 



(1 



Y 



e 



i^X 



AY 
dX 



(117) 



Po cr 
1+0- 



Y- 



AaX 



^X — 
dX 



(118) 

Let us, first, analyze these functions in terms of the 
variable Y. 



A. Qualitative analysis 

1. The model with < ^ < | and —1 < ct < 

Let us assume that the asymptotic solution of 
()116|) at y — oo differs insignificantly from the so- 
lution to (l89l). and let us consider the function 



1 



m 

Y 



(119) 



in which the parameter i?(^) is unknown, but should 
satisfy the condition 



B{0) = 1 + 4a. 



(120) 



This means that we consider the solution in the form 
([57)1 and restrict ourselves by the term y iii the 
parentheses. With mentioned accuracy the function 
X{Y) (ITTO|) satisfies (ITTC| with ([MH) , when 



5(0 = 1- 



4(cT + 30 
3(1 -ho-) 



(121) 



From the qualitative point of view the behavior of 
the DE with ^T^Oatll^-— oois analogous to that 
of the DE with ^=0 and effective parameter a* equal 
to 



g + 3C 
1-3C 



(122) 



It is interesting that o'*=— 1, when a=—l, i.e., 
these two parameters coincide for arbitrary ^ at the 
phantom-crossing point a'=o'*=— 1. 

For the solution of the type ()119|) the term X^ 
behaves as (1-1-^), when y ^ oo. Thus, the con- 
tributions of the last terms in (|117p and (IllSp (pro- 
portional to the parameter S,) can be neglected at 
y — oo. Thus, the main conclusion that the 
Archimedean-type interaction avoids the Big Rip re- 
mains valid atO<^< |, — l<o'<0. 



2. The model with ^ > and a < — 1 

The key equation for this case can be obtained 
from (|116|1 by the formal replacement X ^ —X and 
Y'{X) -Y'{X) with positively defined X. Nu- 
merical calculations show that in this case the regime 
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y — > oo does not exist. To illustrate qualitatively 
such a behavior one can mention the following: the 
regime with Y oo would require that the lead- 
ing order terms in the equation (jll6p are linked by 
Y'e^ +Y=0, which is in contradiction with the re- 
quirement X > 0. Then two scenaria are available: 
first, with Y — > const, second with Y — )■ —oo. In 
the first scenario there is no Big Rip, and, as a rule, 
we deal with the Pseudo Rip. In the second scenario 
the quantity Y with positive initial value Y{1) > 
would have negative derivative, thus, the equation 
(|117[) would be transformed into 



H = - 



A-kG 



II 



Y+^ 
X 



AiX 



dY 



dX 



(123) 

showing explicitly that the derivative of the Hub- 
ble function is negative. This means that at some 
moment t* (finite or infinite) the Hubble function 
will take zero value providing the scale factor a{t) 
to reach its maximal value Omax- In other words, in 
this scenario the Big Rip can nor appear in contrast 
to the model with iy^,=0. 



3. The model with ^ < and cr > — 1 

When ^ is negative, the procedure similar to the 
one used in the case 0<^<|,— l<o'<0 gives the 
asymptotic formulas analogous to (|100p and (|10ip 



i/(y^oo)^^(3-a). 



H'^{Y oo) 



27rG(3- 



Y . 



(124) 



(125) 



which do not contain the parameter ^ in the leading 
order term. However, now we consider the parame- 
ter l-\-a to be positive, and can conclude that we deal 
with the Little Rip, when — 1 < cr < 3, and the solu- 
tion with finite |a(i)|, when cr > 3. In any case the 
Archimedean-type interaction avoids the Big Rip. 



is satisfied with the solution 



H(a;) = H(l) logx, 



(127) 



if the guiding parameters are coupled by the equality 



3^-1 



(128) 



Since i?*/* is positive, one requires that C > 3 + 
\pov*, i.e., this equality can not be satisfied when 
^=0. When the right-hand side of (I128P is posi- 
tive, the logarithmic solution can appear as a result 
of some "fine tuning": one can satisfy (|128p . e.g., 

by varying the initial temperature of the DM, 

2 

which enters the parameter A^, — ^"j; , the argument 
of the function /*(A*) (see (f85|) ). The asymptotic 
DE energy-density is now 



p{x) = p(l) H logx. 



where the relation 



H'(l) 



[p(i)-Po + n(i)] 



(129) 
(130) 



(131) 



is used. The sum ^{x)^ p{x) remains constant 

p{x) + H(a;) - p{l) + n(l) - - — 

The DM energy density E{x) is now constant, i.e., 

E{x) = EJ, , (132) 

and the Hubble function H{x) can be found from 
the equation 

" ^ Ui) 



H\x) 



E,L 



4 , 

— logx 



(133) 



The corresponding scale factor a{t) can be written 
in the form 



a{t) = a* exp 



3i^* 



(t-ty 



(134) 



where the parameters with asterisks are defined as 
follows 



a(to)exp{-^[p(l)+£;(o)]} , 



4-. The special model a——l 

As in the general case described in , asymp- 

totic model admits the solution of the anti-Gaussian 
type, when cr=— 1. Indeed, the corresponding key 
equation 



£,x^n"+{4^-l)xU'+3po 



E,vJ, 



(126) 



t* = to 



I'* 



(o)J 



(135) 



327rG" 

The acceleration parameter —q{t) for the anti- 
Gaussian expansion is positive and exceeds the 
unity: 

-gW^-|. = l+ . >!• (136) 



167rG(t-t*) 



Glearly, this solution illustrates one of the versions 
of the Little Rip. 
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5. Numerical analysis 



In order to confirm the qualitative conclusions 
made for the model with ^ ^ 0, j^* 7^ 0, we ana- 
lyzed numerically the models for the parameters f 
and a belonging to the domains I, II, III, IV displayed 
on Fig.l. These calculations are illustrated by Fig. 2 
and Fig. 3. These figures contain four panels: the 
first one displays the Hubble function H{x), the sec- 
ond panel demonstrates the behavior of the scale 
factor a{t), the third panel describes the rate of the 
growth of the Hubble function H, the fourth panel 
displays the function p{x). 






FIG. 2. The illustration to the model with cr=— 0.1, 
5=-0.1, J/*=l (i.e., for the case -1 < cr < 0, ^ < 
and < 0, which refers to the Big Rip scenario at 

v*—0 according to Fig.l). Clearly, when 7^ 0, i.e., the 
Archimedean-type coupling is switched on, the model 
solution is of the Little Rip type with a — )■ 00, — )■ 00, 
H — >■ const, p — >■ 00, n — > —00. 



VI. DISCUSSION 

When the dark matter is effectively coupled to 
the dark energy by the Archimedean-type interac- 
tion, the late-time evolution of the Universe with 
negative DE pressure is protected from the Big Rip 
singularity, and the Little Rip becomes a typical fate 



H0.4 






10 10 10 

X 



FIG. 3. The illustration to the model with a = —1.01, 
^ = 0.7, u,=l (i.e., for the case a < —1, ^ > and 
a-|-3f > 0, which refers to the Big Rip scenario at v^—O 
according to Fig.l). When 7^ the model solution is 
of the Little Rip type with o — >■ 00, H — >■ 00, H — >■ const, 
p — > 00, n — > —00. 



of the Universe. This is the main conclusion of the 
paper. More detailed discussion includes the follow- 
ing three items. 

1. In the presented model of Archimedean- type cou- 
pling between the dark matter and dark energy the 
DE pressure, H, is the key element of modeling. All 
the model solutions can be divided into three classes 
with respect to asymptotic behavior of the state 
function H(i ^ ioo)=Hoo (we consider both cases: 
too=ts (future finite time singularity) and ioo=oo). 
The first class is characterized by finite Hoo ; the sec- 
ond class relates to Hoo=+oo; we deal with the solu- 
tions of the third class, when Hoc =—00. In [H, 
we focused on the solutions of the first class, and 
have shown qualitatively that the solutions of the 
second and third classes exist for some values of the 
effective guiding model parameters ^, a, po and v^,. 
Three parameters ^, a and po describe the equation 
of state of the dark energy (see ((T^)), i'* is the ef- 
fective Archimedean- type coupling constant. This 
paper is devoted to the analysis of the solutions of 
the second and third classes. When Hoc is finite or 
Hoo=-|-oo, the dark matter is asymptotically decou- 
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pled from the dark energy and its energy density be- 
comes negligible in comparison with the DE energy 
density. We indicated this asymptotic situation as 
the DE domination (the coupling constant j/* is not 
equal to zero, but it becomes a hidden parameter of 
the model). In this case only the appropriate choice 
of the constitutive parameters ^, a and po can pro- 
tect the Universe from the Big Rip singularity. To 
summarize the analytical results for the DE domi- 
nation epochs we prepared Fig.l. 

When ^=0, the Big Rip can be realized on the in- 
terval — 1 < (7 < (or in other words < ^1)- 
It can be indicated as the "classical" Big Rip do- 
mam. When < ^ < i this "classical" Big 
Rip zone contracts along the line a to the interval 
— 1 < cr < 3^— 2-v/3j, and disappears at all when 
^ > i. In this sense, the retardation of the DE 
response to the Universe expansion, can avoid the 
Big Rip, if the relaxation parameter ^ exceeds the 
critical value S,=\- In other words, when the DE 
relaxation time T{t) = -j^j^ is bigger than 3^=^, 

where = VkU^ is the Universe expansion scalar, 
the regime of the Big Rip can not be supported. If ^ 
is negative, the interval — 1 < a < again refers to 
the Big Rip scenario of the Universe expansion. The 
choice of the parameter po can not provide the avoid- 
ance of the Big Rip scenario, but it predetermined 
the type of the Big Rip: the power-law type, hyper- 
bolic or trigonometric ones. When cr > 0, there was 
no "classical" Big Rip at ^=0. The same fact can be 
indicated at ^ > 0. However, if the parameter ^ is 
negative, the Big Rip scenario is possible for cr > 0, 
since, instead of retardation, the DE displays the ac- 
celeration of the response to the Universe expansion. 
When cr < — 1 (or, in other words, — 1 < w < 0), 
there was no "classical" Big Rip at ^=0, however, 
the Big Rip becomes possible at cr < —1 for positive 
relaxation time parameter, ^ > 0. 
2. The most interesting case studied in the pa- 
per concerns the model with asymptotic behavior 
n(t too) "> —00, for which the Archimedean- type 
coupling leads to the effective heating of the dark 
matter component of the dark fluid. Due to the force 
proportional to the four-gradient of the DE pressure 
the DM becomes effectively ultrarelativistic and thus 
plays an active role in the energy redistribution pro- 
cesses inside the dark fluid. The reviving of this sec- 
ond player in the late-time scenario of the Universe 
evolution change essentially the character of expan- 
sion: at I'* ^ the Big Rip scenaria happen to be 
avoided and instead of them the Little Rip scenaria 
become typical for the Universe late-time evolution. 
This avoidance is typical both for the cases ^=0 and 

We would like to emphasize that the so-called anti- 



Gaussian type solutions for the scale factor ait) ap- 
pearing as some specific exact solution in |37| , be- 
come typical asymptotic solutions in the case under 
discussion (see, e.g., (jlOSp and (1134^ '). These anti- 
Gaussian t;^)e solutions correspond to the Little Rip 
scenario [23, [s^, since a{t — 00) — 00, H{t — )■ 
00) — T' 00, |n(t — ^ oo)| — > 00 and p{t — > 00) — > cxi. 
In other words, in this model the life-time of the 
Universe is infinite, the Hubble function and scale 
factor tend to infinity without vertical asymptotes. 

Discussing the mechanism of the Big Rip avoid- 
ance in the presence of the Archimedean-type in- 
teraction, we would like to attract the attention to 
the following feature. When the coupling constant 

vanishes, the DE pressure behaves as the power- 
law function of the scale factor. In terms of dimen- 
sionless variables Y and X (see (|57| ) this power-law 
function is predetermined by the first term in ([55)1 . 
When v^, 0, and the DM contribution to the to- 
tal energy-density gains the same order as the DE 
energy-density, the DE pressure behaves according 
to the super-logarithmic law (|103p . which is typi- 
cal for the Little Rip. Super-logarithmic law for Y 
means that the DE pressure grows more slowly than 
in the case of l'^,—0. In other words, since the DE 
transfers energy to the DM due to the Archimedean- 
type interaction, the rate of the DE pressure growth 
decreases, thus protecting the Universe from the Big 
Rip singularity. 

3. The model under discussion displays that there 
arc three critical values of the guiding parameters. 
The first one is the value cr=— 1, which corresponds 
to the well-known phantom-crossing value ■w—^=~l 
for the DE equation of state parameter w. The 
symptom of criticality is that the term (l+a)^^ sys- 
tematically appears in the key equations and expres- 
sions (see, e.g., (PT|) and ([55]) ). The manifestation of 
the criticality is that the solutions with cr — > — 1 dif- 
fer principally from the ones with a = —1, so that 
we considered the last case as the special one. In 
this sense, one can see some analogy between this 
model and the critical behavior of the magnetic and 
electric fields affected by the gravitational wave: we 
obtained principally different solutions for the case, 
when the refractive index tends to one, n — >■ 1, and 
for the case, when this parameter is equal to one 
identically, n = 1, since the term n^~l appeared 
in the denominators of the sought- for functions [42| . 
The second critical value is z^*=0. Again, the param- 
eter v^, appears in the denominators in the functions 
a{t) and H{t) (see, e.g., (|105p . (|106p . as the symp- 
tom of criticality, and the solutions with I'^.—O (we 
mean the solutions of the Big Rip type) differ prin- 
cipally from the solutions with i^* 7^ (the solutions 
of the Little Rip type). The third critical value is 
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connected with the parameter ^. The value ^=0 can 
be considered as the critical one, since crossing the 
vertical line ^=0 on Fig.l, we see the structural rear- 
rangement of the domains corresponding to the Big 
Rip type solutions. In addition, according to (|83p . 
this parameter appears in the denominator, indicat- 
ing that the situation with ^ — > (there are two 
roots of the characteristic equation) differs essen- 
tially from the situation with = (there exists 
only one root). 

The first application of the model was considered 
in [39| : we hope to present new results concerning 
the cosmological applications of the model in the 



next papers. 
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